Consistent holographic description of boost-invariant plasma 
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Prior attempts to construct the gravity dual of boost-invariant flow of Af = 4 supersymmetric 
Yang-Mills gauge theory plasma suffered from apparent curvature singularities in the late time 
expansion. This Letter shows how these problems can be resolved by a different choice of expansion 
parameter. The calculations presented correctly reproduce the plasma energy-momentum tensor 
within the framework of second order viscous hydrodynamics. 



Introduction 



O PACS numbers: 11.25.Tq, 12.38.Mh 
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f*"** , Gauge/string theory duality [l[ has proved to be a valuable tool in describing the properties of strongly coupled 
C ' gauge theory at finite temperature. Most activity has focused on static situations or linearized dynamics. Only 
recently some progress was made in the studies of strongly coupled nonlinear gauge theory dynamics using string 
■ theory methods 0, H H H S S 1, i, [Tot [Til Il2| . Such studies are currently of great practical importance due to 
experimental investigations of strongly coupled QCD plasma at RHIC and soon at the LHC [13|. While the AdS/CFT 
correspondence is well understood only in the supersymmetric setting, it is believed that this case can capture some 
' essential features of real-world QCD above the deconfinement temperature. 

Experimental data suggest that a boost invariant description of the expansion of the fireball seen at RHIC should 
provide a useful model. This is based on Bjorken's observation [l4j that multiplicity spectra when expressed in proper 
time and rapidity variables are approximately independent of rapidity in the mid-rapidity region. The pioneering work 
of Janik and Peschanski Q established the gravity dual of the boost-invariant flow of AT = 4 plasma in the regime of 
large proper time. These authors showed that asymptotic large time behavior of gauge theory plasma, when matter 
is locally well equilibrated, is given by hydrodynamics. The large proper time expansion is equivalent to a gradient 
expansion of hydrodynamics. In a series of follow-up papers [3|, U, |5[ sub-asymptotic corrections were calculated, 
corresponding to various dissipative terms in the hydrodynamical description. These results were of great interest 
to the heavy ion community, because they provided numerical values of strongly coupled gauge theory transport 
coefficients starting from first principles. Obtaining the energy-momentum tensor of boost-invariant flow up to second 
order also helped to establish the correct theory of causal conformal hydrodynamics . 

The approach explored in 0, H| relied entirely on demanding that coefficients in the expansion of curvature invariants 
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in powers of inverse proper time should be regular. This determined the viscosity coefficient in a way consistent 
with results obtained by other means and for the first time provided information about second order viscous 
hydrodynamics for strongly coupled plasmapl- However, this regularity condition turned out to be violated at third 
order in the large proper time expansion [5[7 It was suggested @, @] that the singularities encountered cannot be 
canceled within the supergravity approximation and they indicate cither the need for additional string theory degrees 
of freedom, or a genuine instability. The intention of this Letter is to readdress this issue. 

Recently another framework describing gauge theory plasma hydrodynamics was developed Q, where the gravity 
dual is determined order by order in a gradient expansion starting from a locally boosted black brane geometry. 
This approach utilizes generalized incoming Eddington-Finkelstein coordinates and yields a manifestly regular metric 
S fill - [HJ • I n @ the energy- momentum tensor was found explicitly up to second order in the gradient expansion; 
the calculation makes it is hard to envisage how singularities could arise at third order, even though such a calculation 
seems technically rather demanding. It is natural to ask whether this experience can be used to shed light on the 
results of [f|. The answer is affirmative. A careful analysis of the late-time boost- invariant geometry employingthe 
incoming Eddington-Finkelstein coordinates shows (as discussed in this Letter) that singularities encountered in [a, E| 
are unphysical and arise due to the choice of expansion parameter. The coordinates used here are suitable, since 
the corrections in the gradient expansion remain valid up to the curvature singularity at the origin. Note that while 
the curvature invariants are diffeomorphism-invariant, both their expansion coefficients in powers of proper time and 
the expansion parameter itself are not. This is the reason why a judicious choice of coordinates is relevant for the 
formulation of a consistent late time expansion. 



Boost-invariant flow from the black brane solution 

Bjorken expansion of Af = 4 SYM plasma is a one-dimensional flow with boost invariance along the expansion 
axis and rotational and translational symmetries in the perpendicular plane [lij ]. Proper time r and rapidity y are 
related to the usual lab- frame coordinates by x° = rcoshy, x 1 — rsinhy. Proper time r is invariant under boosts 
along the collision axis, whereas rapidity y is not. Thus boost invariance implies that physical quantities can depend 
only on proper time r, not on rapidity y. 

Following the ideas of |9| boost-invariant perfect fluid flow can be obtained locally from the 5-dimensional boosted 
black brane solution 

ds 2 = -2u^ L Ar - r 2 - — u^dx^dx" + r 2 (rj^ + u M n„) dx^dx", (1) 

where it M is the boost velocity parameter and b is a dilatation parameter related to the black brane temperature 
T by b = 1/ttT. The key ingredient of this approach is the introduction the incoming Eddington-Finkelstein time 
coordinate. At the boundary this coordinate reduces to the usual Minkowski time. Gauge/gravity duality maps 
physics at the boundary to the bulk of the AdS spacetime. For an effective description at long wavelengths (relativistic 
hydrodynamics) the map is provided by decreasing r while keeping the Eddington-Finkelstein time coordinate fixed 
[1 01 ] . For a Bjorken expansion it is natural to use proper time r instead of the usual Minkowski time, so an analogous 
Eddington-Finkelstein type proper time coordinate f is introduced. Specifically, u = d T at the boundary, but is now 
taken as u = df in the bulk of AdS space. Furthermore, for a boost invariant flow the temperature T is asymptotically 
proportional to r -1 / 3 , which translates to b = 3 1 / 4 2 _1 / 2 f 1 / 3 (to ensure agreement with 0]). Thus finally 

ds 2 = -r 2 (l i4r-r)df 2 + 2dfdr + r 2 f 2 dy 2 + r 2 dx\. (2) 

It is straightforward to verify that metric is related to the Janik-Peschanski metric 

d ^P = ^ (- -^r^ dr 2 + r 2 (l + ^ ) d, 2 + (l + ^3 ) d*i + dz 2 1 (3) 
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by a coordinate transformation 
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Because of the nontrivial dependence of f on r the limits r — > oo and f — > oo (corresponding to equilibration of gauge 
theory plasma) differ in the bulk. Moreover the metric @ is regular and invertible up to the black brane singularity, 
which is not the case with ([3|). These observations are crucial for the present approach to work. Note also that the 
relation between r and f is singular when z — 3 1 / 4 t 1 / 3 . This is precisely the locus where the singularities found in 
Q were encountered. 

The metric ([2]) is not an exact solution of Einstein equations Rm n + 4G_a/ tj = - there are subleading corrections 
coming from derivatives of the velocity u and temperature T. They correspond to the gradient expansion of the 
boundary energy momentum tensor [17] . For a boost-invariant flow the energy- momentum tensor < T^ v > is deter- 
mined by the energy density e(r) Q. In the gradient expansion each covariant derivative V« is damped by jKp (L 
is the characteristic length scale of a perturbation and T is the fluid temperature). Because boost-invariant flow is 
characterized by T ~ r -1 / 3 + . . . and moreover Vw ~ r — 1 , we expect the following expansion of the energy density 
<(t) M 

<t) = ^7s I 1 - **>^ + {^nl ^or° + |a?) ^ + . ••} , (5) 

where 770 and r^} are the viscosity and relaxation time coefficients, whereas Aj is a new transport coefficient introduced 
in (l5j . The above expansion is written in terms of r, since the energy density is defined at the boundary. 

The apparent curvature singularities in AdS encountered in [j| appear at third order in the large r (gradient) expansion. 
It is difficult to check what happens for a general flow at this order. However the situation is much simpler for a 
boost-invariant flow, since all the symmetries can be imposed from the outset. This leads to the following ansatz for 
the metric [2l| 

ds 2 = G MN dx M dx N = -r 2 N{f, r)df 2 + 2dfdr + r 2 r V (f ' r) dy 2 + r 2 e c(f - r) dx 2 ± . (6) 

Introducing the scaling variable v — r ■ f 1 / 3 in analogy with what is done in 0, d, & H, B 0, II] one obtains the natural 
expansion of the metric components in f~ 2 / 3 on the gravity side 

N(f,r) = A(v)exp(j2 a ^ f ~ 2k/3 )> 

\A;>0 / 
e b(f,r) = B ( w ) exp I J2b k {v)f- 2k / 3 ), 



\fc>0 



C( V )exp [J2 c k(v)T- 2k ^ . (7) 



\k>0 



Note that the scaling variable v introduced here is different from the one used in 0, H, H]> H, H, 0; Q , and so the scaling 
limit considered here must be regarded as different. To obtain a uniform expansion of the Einstein equations Emn = 
Rain + ^Gmn = one needs to rescale them (see according to E = (f 2 / 3 Eff, E? r , j^E rri ^j^E yy , t 2 / 3 E x±x± ) . 
This leads to 

E(f,r) = E (r ■ f 1 ' 3 ) + -^Mr ■ r 1 ' 3 ) + ~^E 2 {r ■ r 1 ' 3 ) + ^E 3 (r ■ r 1 ' 3 ) + ... (8) 

The curvature invariants (e.g. TZmnopT^- MNOP ) defined recursively in [f| can be likewise expanded. The crucial 
difference between the present approach and the one introduced in is that the expansion parameter involves f 
instead of r. Einstein equations can be solved order by order in f~ 2 l 3 expansion starting from 

Mv) = 1- A, 

3v 4 

B(v) = C(v) = l, (9) 

which simply reproduces the boosted black brane solution ©. Thus the zeroth order solution entails large but finite 
f ; it is to be expected that the singularity at r — should be shielded by an event horizon, however it is difficult to 
demonstrate this explicitly^. 
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Gravity dual of the gradient expansion 

The equations of motion (|SJ) at a given order k are a system of ordinary second order differential equations for the 
3 functions (v), bk {v) and (v). Each solution involves two integration constants. On the other hand, two of 
the equations of motion are constraints. At each order k > one of the constraints fixes one of the integration 
constants appearing at that order, and the other one fixes an integration constant left undetermined at order k — 1. 
The 4 remaining integration constants can be fixed order by order by imposing metric regularity (up to the usual 
black brane singularity at v = [9j). It turns out that the potential singularity is located only at v = \/2/2> 1 l i ; 
thus the functions bk (v) , Cfc (v) must remain finite as v — > V2/3 1 / 4 . In case of a,k (v) the requirement should be 
that the product with A(v) must be finite. However one can take advantage of residual diffeomorphism invariance 
[HI preserved by the ansatz (jSJ), which can effectively be fixed by requiring that a,k (v) itself be regular. Furthermore 
asymptotic AdS behavior of the metric requires that these functions vanish as v — > oo (in the late proper time regime) . 
These conditions together with the constraints fix 5 of the 6 integration constants at a given order k > and lead to 
a regular metric with no poles or logarithmic singularities apart from v = 0. As an example, the first order solution 
(dual to viscous hydrodynamics [3j) reads 

_ 2 2-3- 1 / 2 + 2 1 / 2 3- 1 / 4 v + v 2 
ai 3 ' (2V23-1/4 + „) . (2 . 3-Va + v 2 ) ' 

n V2 /3 1 / 4 \2V2, V2 f V2 \ 1 , / 2 2 \ 

Oi = —f= — ; 777 arctan — =-« ^-rr logv H — -77 log —777 + v M t= — t log — ttt + u ■ (10) 

V23 3 / 4 3 3 / 4 VV2 / 3 3 / 4 B 3 3 / 4 B \ y 3 1 / 4 y x/23 3 / 4 \3 1/2 / 

and ci = —bi/2. Higher order formulae (up to the third order) can be found in a Mathematica notebook available 
online |l9j |. Holographic renormalization correctly reproduces the energy density for the boost-invariant flow up to 
second order in derivatives [j| [H[ . 

Absence of singularities and relation to Fefferman-Graham coordinates 

The assumption of non-singularity of coefficients of curvature invariants in the late proper-time expansion was crucial 
in establishing transport properties of J\f = 4 SYM plasma in The present approach to boost- invariant flow starts 

from a manifestly regular metric in the leading order (no logarithmic and power-like sing ularities at v = V2/3 1/4 ) 
and produces regular solutions up to the third order. Since the components of the metric as well as its inverse are 
regular (as well as their derivatives), all curvature invariants are non-singular. Indeed, from ([6]) it follows that the 
non-vanishing components of the inverse are G rr = r 2 e a ( T > r \ G rT = 1, G yv = r~ 2 f~ 2 e~ b ( 7 '' r \ G ±J - = r~ 2 e~ c< -- T > r \ If 
e -a(r,r) na( j Deen present, singularities would have appeared as a consequence of ©• 

It is natural to ask how these results are related to those obtained using the original approach of [2j, |5( . Clearly, they 
should be related by a coordinate transformation order by order in the large proper-time expansion: 




It is straightforward (though tedious) to determine this transformation explicitly; the results up to third order are 
online [19(. This transformation defines a map between the two approaches and explains how the results obtained in 
H were correct despite the singularities which were encountered there. 

Summary 

This Letter shows how to construct a consistent gravity dual to boost-invariant flow of Af — 4 SYM plasma in 
a gradient expansion. This solution reproduces known results up to second order viscous hydrodynamics A 
manifestly regular metric up to third order was described and arguments were given why all curvature invariants have 
non-singular expansions in f ~ 2 / 3 . The relation with the approach based on Fefferman-Graham coordinates was made 
explicit by the coordinate transformation discussed above. 

In conclusion, the AdS/CFT dual to boost-invariant flow can be realized within the supergravity approximation and 
is completely free of singularities apart from the black brane singularity at r = 0. 
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Note added: Shortly after the original version of this Letter was posted to arXiv.org a very interesting paper 
(Tl | appeared which deals with the same subject and along similar lines as the approach disussed here. Apart from 
demonstrating the presence of an apparent horizon, the authors also give an all-orders argument for the absence of 
singularities in the large proper time expansion. 
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